Factorial cumulants reveal interactions in counting statistics 
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Full counting statistics concerns the stochastic transport of electrons in mesoscopic structures. 
Recently it has been shown that the charge transport statistics for noninteracting electrons in a 
two-terminal system is always generalized binomial; it can be decomposed into independent single- 
particle events, and the zeros of the generating function are real and negative. Here we investigate 
how the zeros of the generating function move into the complex plane due to interactions and 
demonstrate that the positions of the zeros can be detected using high-order factorial cumulants. 
As an illustrative example we consider electron transport through a Coulomb blockade quantum 
dot for which we show that the interactions on the quantum dot are clearly visible in the high-order 
factorial cumulants. Our findings are important for understanding the influence of interactions on 
counting statistics, and the characterization in terms of zeros of the generating function provides us 
with a simple interpretation of recent experiments, where high-order statistics have been measured. 

PACS numbers: 02.50.Ey, 72.70.-|-m, 73.23.Hk 
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I. INTRODUCTION 

Full counting statistics (FCS) has-Jpcen a topic of ac- 
tive research for nearly two decades. ElH FCS describes the 
statistics of charge transport through mesoscopic conduc- 
tors [Fig. 0(a)] and is expected to provide more informa- 
tion about the physical processes inside a conductor com- 
pared to what is available from the mean current and the 
shot noise only. Within this framework it is natural to ask 
how such additional information can be extracted from 
the high-order statistics and what quantities or measures 
are most suitable to this end. This core problem of FCS 
constitutes the main focus of the present work. 

On the experimental side, FCS has recently gained con- 
siderable impetus due to a number of measurements of 
high-order statistics in nanoscale systems. While earlier 
experiments were restricted to the first few moments or 
cumulants of the current, high-order cumulants of the 
charge transport are now becoming experimentally ac- 
cessible. The fourth and fifth cumulants of the current 
have recently [been measured both in Caulomb blockade 
quantum dotsO and in avalanche diodesJj Quantum dots, 
in particular, have emerged as useful sources of high- 
accuracy counting statistics. In these systems, the low 
(kilo-hertz) tunneling rates allow for detection of individ- 
ual electrons in real-time using a nearby quantum point 
contact whose conduC|tancc is sensitive to the charge oc- 
cupations of the dots.fflia Although quantum effects are 
typically suppressed at the long time scales characteriz- 
ing the charge transport, quantum dots provide a unique 
setting to experimentally test theoretical predictions for 
high-order statistics. Remarkably, time-dependent cu- 
mulants of the transferred charge beyond the 15th or- 
der have recently been measured in single-electixin. trans- 
port through a Coulomb blockade quantum dotLTlj and a 
wealth of statistical data is now available. 

On the theory side, several techniques have been devel- 
oped for calculating the statistics of transfesred charges. 
In the seminal work by Levitov and LesovikJ^I the FCS of 
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FIG. 1. Nanoscale conductors, (a) Generic transport setup 
consisting of a nanoscale conductor (gray) connected to source 
and drain electrodes. Electron flow in both directions (ar- 
rows) is allowed. The probability distribution for the number 
of electrons n that have been collected in the drain electrode 
during the time span [0,t] is denoted P{n,t). (b) Single- 
level quantum dot coupled via tunnel barriers to source (S) 
and drain (D) electrodes, kept at temperatures Ts and Tn, 
respectively, and chemical potentials us and fio- The single- 
level energy of the quantu m dot is Sd and the on-site Coulomb 
interaction is U (see Sec. Ill A). 



noninteracting electrons propagating coherently through 
a conductor was expressed by a determinant fosp 
containing the scattering matrix of the problem.Elll2r 
In many-channel conductors the statistics is predomi- 
nantly classical and can b^ described using a Langevin- 
Boltzmann equation.E3Ej A powerful and elegant formu- 
lation in this scmiclassical regime is the stochastic path 
integral approach, which was introduced for FCS in the 
pioneering works by Pilgram et aZ.EjtZl For interacting 
systems, the FCS can be related to a generalized master 
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equation describing the charge tr|aiisBortt3i22l or obtained 
using Keldysh Greens functions Jlj'EHl 

The central object in the theory of FCS is the generat- 
ing function (GF) for the probabihty distribution t) 
of the number of transferred electrons n. Its analytic 
properties as well as its symmetries as a function of volt- 
age, temperature, and counting ficlda_which, for exam- 
ple, lead to fluctuation relations, E3c3 are therefore of 
fundamental interest. If the transport process consists 
of independent, elementary events, the GF may be fac- 
torized according to these. This has been pointed out by 
Vanevic, Nazarov and Belzig, who identified the elemen- 
tary transport processes in a qiiaxitum conductor driven 
by a time-dependent voltage.EZiEEl Importantly, Abanov 
and Ivanov have shown recently that for noninteracting 
electrons in a two-terminal setup, the GF can always be 
factorized into single-particle events and the zeps-of the 
GF correspondingly lie on the negative real axis.ESE^ This 
form of the counting statistics has been dubbed general- 
ized binomial by Hassler et alS3 Interactions, however, 
can change these properties and cause the zeros to move 
into the complex plane.E2l As evident from this discus- 
sion, the zeros of the GF are of crucial importance for 
understanding the physical mechanisms that determine 
the transport statistics. However, given a complicated 
GF or even actual experimental data, it is not clear how 
one can tell, whether or not the statistics are generalized 
binomial and if the charges interacted inside the conduc- 
tor. 

Traditionally within FCS, the probability distribution 
P{n, t) has been characterized by its cumulants, from 
which one could hope to extract infornijition about the 
physical processes inside the conductor.cl The cumulants 
are related to derivatives of the GF. Several theoretical 
studies, however, have found that the high-order cumu- 
lants ten d to nscillate strongly as some system parameter 
is varied.c3E3^Zl Recently, this phenomenon has been ex- 
plained on general grounds and it has been shown that 
the high-order cumulants for almost any non-Gaussian 
distribution should|X)scillate as functions of basically any 
system parameter,Ll following wei-k by Berry on high 
derivatives of smooth functions.E§ The fact that these 
oscillations are so generic could indicate that the high- 
order cumulants are fragile objects from which it is dif- 
ficult to extract information about a particular system: 
The information in the high-order cumulants seems to 
be masked by the oscillatory behavior that should be 
present in almost any system. However, as we explain 
below, such oscillations do not prevent high-order statis- 
tics from bearing information. 

It is the purpose of this paper to argue that factorial 
cumulants provide an alternative characterization of the 
probability distribution P{n, t) that is particularly useful 
for describing the high-order statistics of nanoscale trans- 
port. In mesoscopic physics, factorial cumulants have so 
far received only limited attention, except in stijdies of 
photons emitted from a quantum point contact,E2l and, 
very recently, in connection with current fluctuations and 



entanglement entropy.0 However, as we show, the high- 
order factorial cumulants provide a simple description of 
generalized binomial statistics and they directly reflect 
the zeros of the GF. In fact, for noninteracting systems, 
the high-order factorial cumulants do not oscillate (un- 
like the ordinary cumulants), no matter what parameter 
is varied. In contrast, if a high-order factorial cumulant 
oscillates as a function of some parameter, it cannot be 
describing generalized binomial statistics and the elec- 
trons must have interacted inside the conductor. We il- 
lustrate these points with a model of transport through 
a quantum dot [Fig. |l|b], for which we show how inter- 
actions on the dot cause the zeros of the GF to move 
into the complex plane as clearly seen in the factorial 
cumulants. 

The paper is organized as follows: In Sec. ^ we in- 
troduce the general framework of FCS, including the GF 
and the moments and cumulants, as well as their facto- 
rial counterparts. We show that the factorial cumulants 
are particularly simple for generalized binomial statistics, 
and we give an asymptotic expression for the high-order 
factorial cumulants, which is directly related to the zeros 
of the GF. We show that the behavior of the high-order 
factorial cumulants changes drastically as the zeros move 
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into the complex plane due to interactions. In Sec 
we illustrate our findings using a model of electron trans- 
port through a quantum dot, which is weakly coupled to 
source and drain electrodes. At low bias voltages, the 
quantum dot can only be empty or singly occupied. Th: 
corresponds to the recent experiment by Fricke et al. 
and we show how the measured oscillations of the 15th 
cumulant as a function of time and coupling to the leads 
can be explained by the motion of the zeros of the GF. 
In contrast, the factorial cumulants do not oscillate. As 
the voltage difference between the leads is increased, the 
quantum dot can be occupied also by two electrons at a 
time, and the on-site Coulomb interaction now strongly 
affects the FCS. The zeros of the GF move into the com- 
plex plane, which is clearly visible in the high-order fac- 
torial cumulants. Finally, Sec. ^ contains our conclud- 
ing remarks as well as a number of open questions and 
directions for future research. Appendixes A| and re- 
spectively, describe methods for calculating high-order 
(factorial) moments and cumulants at finite times from 
a master equation and for extracting the position of the 
zeros from the high-order factorial cumulants. 



II. FULL COUNTING STATISTICS 

We consider a generic transport setup in which a 
nanoscale conductor is connected to source and drain 
electrodes [Fig. 0a]. A bias voltage (possibly time de- 
pendent) between the electrodes drives electrons through 
the conductor, and we denote by P{n,t) the probabil- 
ity distribution of the number of electrons n that have 
traversed the nanoscale conductor during the time span 
[0,t\. Depending on the direction of the tunneling events 
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between the source and the drain electrodes, the number 
of transferred electrons n can be either positive or neg- 
ative. For concreteness we define electron flow into the 
drain electrode as the positive direction of the current. 
Examples of nanoscale conductors used in recent count-, 
ing statistics experiment|S-|include tunnel junctionsJid'LHI 
quantum Doint contacts,E3 double quantum dots,0 and 
nanowires,c3 but for the discussion in this section it is 
not necessary to specify the details of the system. 

A. Generating function 

The object of main concern in FCS is the GF, defined 

as 

n 

The GF encodes the full information about the prob- 
abilities P(n, t) and it allows us to introduce (facto- 
rial) moments and (factorial) cumulants in the follow- 
ing. As mentioned in Sec. |, calculations of the GF 
can be approached with several different techniques, 
depending on the specific system at hand. Interest- 
ingly, however, Abanov and Ivanov have recently shown 
that the counting statistics for a two-terminal scatter- 
ing problem involving noninteracting electrons can al- 
ways be decomposed into independent single-particle 
events for as m i fo pn of the scattering matrix and at any 
temperature£al23 Formally, this caja.|-he expressed as a 
factorization of the GF of the formESEJ 

^(^,oT="T^"''n^'(-'*)' (2) 

binomial 

wlicM' Qi{z,t) = 1 — Pi + piZ is the binomial GF corre- 
sponding to a single-particle event occurring with prob- 
ability < Pi < 1, depending on time t as well as all 
other parameters of the system. The factor z~'^ cor- 
responds to a deterministic background charge transfer 
Q = J2iPi ~ (■"■) — opposite to the positive direc- 
tion, where (n) is the mean value of the total transferred 
charge. For unidirectional transport, we have Q = 0. 

A statistical distribution given by Ea- (H) has been 
dubbed generalized binomial statistics O Importantly, 
the result of Abanov and Ivanov implies that if a GF 
cannot be factorized as in Eq. (^) with real probabili- 
ties Pi, it cannot be describing noninteracting electrons. 
We note, however, that the opposite is not true: even 
in the presence of intetactions, the statistics may still be 
generalized binomial .EJ 

B. Moments and cumulants 

We now turn to the moments and cumulants of P{n, t), 
which arc commonly used to characterize the probability 
distribution. The moments (n™) can be found from the 



moment generating function which is obtained from the 
GF [Eq. (|^)] via the substitution z e^, 

A^(z,^) = a(e^^)=^P(r^,^)e"^ (3) 

n 

The moments of the transferred charge are given by the 
derivatives of the moment generating function with re- 
spect to z, evaluated at z = 0, 

(n'")(i) =^n"P(n,t) = dTM{z,t)\^_^,. (4) 

n 

The cumulant generating function (CGF) is defined as 

5(z,t) = ln[X(z,i)] =ln[g(e^^)] (5) 

and the cumulants are similarly defined as derivatives of 
the CGF at z = 0: 

{{n'^)){t)=dTS{z,t)\^^,. (6) 

The first cumulant ((ri)) = (n) is the mean of n, the 
second ((n^)) = (n^) — (n)^ is the variance, and the third 
({n^)) = ((n— (71))'^) is the skewness. For a Poisson distri- 
bution all cumulants are equal to the mean {{n"^)) = (n), 
while only the first and second cumulants are nonzero for 
a Gauss distribution, that is, {{n"'')) = for to > 2. 

C. Factorial moments and factorial cumulants 

A complementary characterization of the probability 
distribution is provided by the facto rial j noments and 
the corresponding factorial cumulants.EJcZl The factorial 
moments are defined as follows 

= (n(n-l)---(n-m + l)). (7) 

It is easy to show that they are generated by the function 

Mriz, t) = g{z + 1, i) = ^ t){z + 1)", (8) 

n 

obtained from the GF [Eq. (|])] via the substitution z — !> 
z + 1. Thus 

{n-^)Fify = d^MF{z,tX^^, (9) 

and similarly to the cumulants, the factorial cumulant 
generating function (FCGF) is defined as 

5i.(z,t) =ln[A^j^(z,i)] =lii[S(;2 + l,i)], (10) 
whose derivatives at z = deliver the factorial cumulants 

inn)F(t)^K?SF{zM.^^- (11) 

The first two factorial cumulants are ((n)) p = (n) and 
{{ri^))F — {n^) — (n)"^ — (n). For a Poisson distribution. 
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only the first factorial cumulant is nonzero and {{n^))F = 
for 771 > 1. 

The factorial cumulants can be expressed in terms of 
the ordinary cumulants via the relations 

{{n'))F = in')) in')). (12) 
in^))F^in'))-'iin'))+2in')), 

which, for arbitrary order m, readS 



(13) 



where s (rn,j) are the Stirling numbers of the first kind. 
They can be generated fmjn the relation [In (1 + x)y = 

In the case of generalized binomial statistics [Eq. (H)] 
the expression for the factorial cumulants becomes par- 
ticularly simple: 



((n-))f^«'="°="'='(-l)™-i(m-l)! 

binomial 



(14) 



This expression provides us with a direct test of whether 
or not a statistical distribution can be factorized into 
independent single-particle events as described by Eq. 
(^. In particular, for unidirectional transport (where 
Q = 0), the factorial cumulants must have alternat- 
ing signs as functions of the cumulant order m due to 
the factor (—1)'""-'^, if the statistics is generalized bino- 
mial. We remark that it is straightforward to extend this 
analysis to bi-directional transport. In that case, the 
quantity ((n"))^ - (-l)'""i(m - l)!(n) = (-l)™"i(m - 
1)! (p™ — Pi) has alternating sign as a function of m 
for generalized binomial statistics. In this work we use 
Eq. (Q) to test whether the counting statistics of charge 
transport through a nanoscale conductor is generalized 
binomial. Importantly, factorial cumulants are measur- 
able and the test is consequently immediately applicable 
to experimental data. 



D. High-order (factorial) cumulants 

As the zeros of the GF move into the complex plane 
due to interactions, the behavior of the high-order facto- 
rial cumulants changes compared to that of generalized 
binomial statistics [Eq. (Q)]. To understand this, we 
note that both the CGF and the FCGF, defined in Eqs. 
(^ and (|lo|), respectively, have logarithmic singularities 
corresponding to the zeros of the GF. These singularities 
determine the high-order asymptoticg|p£|thp ordinary and 
the factorial cumulants, respectively.LrEilHca In this work 
we consider cases for which the (F)CGF has only loga- 
rithmic singularities, as is typical at finite times, but in 
the long-time limit it may have branch-point singulari- 
ties, as we discuss at the end of this section. 



Wc first note that the CGF (or the FCGF) close to a 
logarithmic singularity Zj with degeneracy aj behaves as 



S(^p-f{z,t) ~ aj ln(2j — 2), z close to zj, (15) 
with corresponding derivatives reading 
(m - 1)1 

9"5(i=^)(z,i) ~ -Qfj- — — — -, z close to Zj. (16) 

(Zj z) 

The high-order derivatives evaluated at 2; = 0, that is, 
the high-order (factorial) cumulants, can then be approx- 
imated as a sum over all singularities 



in"^))^P)=dTS^F){z,t)U^o 
^-(m-1)!^, 



according to the first Darboux approximation@@ 

Equation (|l7| ) shows that the high-order (factorial) cu- 
mulants are determined by the singularities closest to 
z = 0, which dominate the sum for large m. Relative 
contributions from other terms are suppressed with the 
relative distance to z = and the power mB If the closest 
(nondegenerate) singularity, denoted zq, lies on the neg- 
ative real axis such that zq = |zo|e"^, Eq. ( p^ ) becomes 
particularly simple for large m and reduces to 



){F) 



(-l)(™-i)(m-l)!/|zor 



(18) 



This is the case for high-order factorial cumulants cor- 
responding to generalized binomial statistics [Eq. (^]. 
For unidirectional transport {Q = 0) the corresponding 
FCGF is 

SFiz,t)'-=''''''Y.^n{l+p,z), (19) 

binomial ^ — ^ 



which has logarithmic singularities at Zj = —1/pj < — 1. 
The high-order factorial cumulants arc then dominated 
by the singularity Zmax = ^1/Pmax closest to z = 0, 
where Pmax is the largest probability among the pi's. We 
thus find 



in"^))F ^^^4"^'='^ (-1)- 

binomial 



(20) 



for large m. For the above example, this conclusion could 
also have been reached directly from Eq. ( p^ (for Q = 0). 

The singularities, however, do not always lie on the 
negative real axis, but in general they come in complex- 
conjugate pairs, ensuring that the (factorial) cumulants 
arc real. In case only a single complex-conjugate pair of 
logarithmic singularities, zq and Zg, is closest to z = 0, 
Eq. ( p7| ) simplifies to 



lnn)iF) ^ 



2(m- 1)! 

i 



COS (marg [zq]) (21) 
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for large m. In the case where arg [zq] = tt, this reduces 
to the right-hand side of Eq. (|l^) muhiphed by 2, since 
the two singularities are then degenerate. 

Interestingly, our analysis shows that high-order fac- 
torial cumulants corresponding to generalized binomial 
statistics have a sign that is determined solely by the 
order m via the factor In contrast, if the 

statistics is not generalized binomial, the factorial cumu- 
lant of a given order m will oscillate as a function of 
any parameter that changes the position of the singular- 
ities due to the factor cos (to arg [zq]), which also causes 
trigonometric oscillations as function of the order to. 

In this section, we have analyzed the situation where 
the (F)CGF has logarithmic singularities due to zeros of 
the GF. It is, however, well known that the (F)CGF can 
have, for example, branch-point singularities in the long- 
time limit, and a generalization of the above analysis is 
necessary to treat such cases. The overall conclusions, 
however, remain intact also for branch-point singulari- 
ties, and since we mainly consider finite times in this 
work, we do not encounter such situations. Instead, we 
refer the interested reader to Refs. |^, |3^, and ^ and, in 
particular, section IV of Rcf. ^ for a more general anal- 
ysis of high-order derivatives and (factorial) cumulants 
for (F)CGFs with branch-point singularities. 

In the next section, we illustrate how the statistics of 
charge transport through a quantum dot due to interac- 
tions can change from being generalized binomial, with 
factorial cumulants given by Eqs. (^4|), to a different sta- 
tistical distribution, with high-order factorial cumulants 
governed by Eq. (pT|). 

III. COULOMB BLOCKADE QUANTUM DOT 
A. Model 

We consider electron transport through a quantum 
dot (QD) with a single spin-dcgencratc level coupled to 
source and drain electrodes [Fig. |l|b] . The energy of the 
level is denoted by Sd and U is the on-site Coulomb in- 
teraction. The Hamiltonian for the coupled system reads 

H = Hd + HT + Hr, (22) 

where 

Hd^ ed{n^ + hi) + Uh^n^ (23) 

is the Hamiltonian of the QD, tunneling between the QD 
and the leads is given by the term 

Ht= {to.kcl^X + h. c.) (24) 

a=S,D 

and the source (a = S) and drain (a = D) electrodes are 
described as reservoirs of free electrons with energy Eaka 



Hr= ^ikJaka- (25) 

k,(T. 
a=S,D 

Here we have defined the fermionic operators dj. (do-), 
which create (annihilate) electrons with spin a on the 
QD, and the corresponding spin-resolved occupation 
number operators are = d^da, u =t,]r. The oper- 
ators c^j,^ and Caka create and annihilate, respectively, 
electrons with momentum k, spin cr, and energy Eaku in 
the source {a = S) or drain (a — D) electrodes. We 
assume that the tunneling matrix elements tak are in- 
dependent of spin and consider the situation without an 
applied magnetic field. Both of these assumptions can 
be lifted, although such extensions of the model are not 
considered here. 

In the following, we consider weak coupling between 
the QD and the leads. Electron transport through the 
QD can then be described by a master equation for tran- 
sitions between different many-body eigenstates of the 
QD. The eigenstates of the isolated QD corresponding 
to Eq. (H) arc |0), | t), I i), and |2), where the first 
and the last cigenstate correspond to the QD being oc- 
cupied by zero or two electrons, respectively, while the 
other two correspond to the QD being occupied by a 
single electron with spin a =t, i, |£fis«ectively. We now 
define the n-resolved probabilitiesllSo po{n,t), p^{n,t), 
Pl{n,t), and p2{n,t) for each of the eigenstates to be 
occupied, while n electrons have been collected in the 
drain during the time span [0,t\. Since tunneling is spin 
independent, it is equally probable to occupy each single- 
electron spin state, Pi^n, t) = (n, t), and we can define 
Pi{n,t) = pi(n,t) +pf{n,t) and collect the probabilities 
in the vector 

\p{n,t))) = [po{n,t),pi{n,t),p2{n,t)f, (26) 

where we use double-brackets to avoid confusion with the 
quantum states of the Hamiltonian. We also define ((0| = 
[1, 1, 1], allowing us to express the probability P{n,t) as 
the inner product P(n,<) = {{Q\p{n,t))) . The GF is then 

g{z,t)^mg{z,t))), (27) 

where \g{z,t))) = T,^z'^\p{n,t))). We note that 
\g{l,t))) = [po{t),pi{t),p2{t)]'^ contains the probabili- 
ties Pi of occupying the QD with i = 0,1,2 electrons 
independently of the number of transferred electrons n. 
As a consequence of probability conservation, we have 

gil,t) = ((0|g(l,0)) =E„^Ki) = 1 for all t. 

The time dependence of the GF is determined by the 
dynamics of \g{z,t))). We find the time evolution of 
\g{z,t))) by setting up a master equation for the tran- 
sitions between the eigenstates of the QD. In the weak 
coupling regime, the transition rates can be found us- 
ing Fermi's Golden rule, treating the tunneling Hamilto- 
nian in Eq. (^) as the perturbationrj Working in the 
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wide-band limit and assuming a constant tunneling den- 
sity of states in both leads, we define the bare energy- 
independent tunneling rates Fq = 2Tr\tka\'^'Da, where Va 
is the density of states in lead a ~ S,D. The master 



equation for \g{z,t))) then reads 

dMz,t)))^M{zMz,t))), (28) 
with the z-dependent rate matrix 



M(z) 



-2(r54°) + rong)) 
2(rs4'' 





zTDil-n^D^)+rsil-n's') 



(0)n 







-roil 



(U) 

z^^Vun 



(U) 
D 



^f) 2{Ts{l~n^s>) + zTn{l-n';^>)} 



-2{rs(i 



AU) 



) + ri5(i 



)} 



Here, we have introduced the Fermi functions of the 
leads, kept at electron temperature Tq, and chemical po- 
tential fia , evaluated at the energy Sd + E 



AE) 



1 



S,D. 



(29) 



In the rate matrix above, only na"^ and 
which we use to parametrize the applied voltage biases 
and the temperatures of the electrodes in the following 
subsection. We note that the off-diagonal elements of 
'M.{z) include factors of z and z^^ multiplying the rates 
corresponding to processes that increase or decrease, re- 
spectively, by 1 the number of electrons that have been 
collected in the drain.E3 

We find the GF by formally solving Eq. (^) for 
\g(z,t))). To this end, we need to define an appropri- 
ate initial condition. We assume that the occupations 
of the QD eigenstates have reached their steady state 
at t = when counting begins and P{n,t = 0) = Sn.o- 
We then have \g{z,t = 0))) = |0)), where the station- 
ary state |0)) of the QD is given by the unique solution 
to M(z = 1)|0)) = 0. Wc thereby obtain the following 
compact expression for the GF: 



(u) 

Ha appear, 



,M(z)t| 



(30) 



This is a general and formally exact result, but in prac- 
tice, given a rate matrix M(z), it may not be possible to 
obtain a simple, closed-form expression for the GF. Fur- 
ther complications arise when trying to calculate (facto- 
rial) moments and (factorial) cumulants, since derivatives 
of the GF with respect to z are required. To calculate 
high-order (factorial) cumulants at finite times we have 
thus developed the method described in Appendix ^. 



B. Results 

We concentrate on the situation, where the leads are 
voltage biased such that the energy level of the QD is 
well above the chemical potential of the drain, I^d 
or n 



(0) 
D 



0, and electrons cannot tunnel back 



is below the chemical potential of the source electrode, 
£d ^ Ms or n^g ^ ~ 1. Under these voltage conditions, 
electron transport takes place from source to drain via 
the QD. The singly occupied state of the QD always par- 
ticipates in transport, but the doubly occupied state only 
becomes populated for nonzero values of n^^\ In the 
following we study the charge transport statistics as a 



function of 



AU) 



that is, we vary the bias in the source 



electrode such that 



AU) 



takes values between and 1. 



Low bias 



We first analyze the low-bias regime n^g^^ ~ 0. This 
corresponds to the situation recently investigated in the 
experiments described in Rcfs. In Rcf. ^it was found 
experimentally that the high-order cumulants oscillated 
as functions of the dimcnsionlcss time 



and the asymmetry 



2rst 



2Ts - Td 



2r, 



T 



D 



(31) 



(32) 



Here, we have included a factor of 2, corresponding to the 
two spin species in our model. In Fig. ^ we show nu- 
merically exact results for the high-order cumulant {{n^^)) 
corresponding to the measurements in Rcf. ^ Our calcu- 
lations agree well with the experiment and reproduce the 
clear oscillations as functions of the dimensionless time r 
and the asymmetry a. To understand the experimental 
and numerical results we next consider the zeros of the 
GF. 

At low voltages, the doubly occupied state of the QD 
remains unpopulated in the stationary state and wc can 
set p2{n,t) ~ 0. This implies a significant simplification 
of Eq. (pq) , which in this case can be written as 



(33) 



d 


9o{z,t) 




■-2rs zTd' 




go{z,t) 




dt 






2Ts -Td 




_9i{z,t) 




where Qi 


{z,t) = 




z'^Piin^t), i -- 




0, 1. We can 



into the QD from the drain. At the same time, the level explicitly evaluate the GF and we obtain the analytic 
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a 



FIG. 2. (Color online) High-order cumulants in the low-bias regime, ~ 0. (a) The cumulant {{n^'')) as a function of 

the asymmetry a [Eq. (^2|)] and the dimensionless time r = 2Vst. The high-order cumulant oscillates as a function of both 
parameters, in agreement with the experimental results of Fricke et al. in Ref. |^ (Note: We do not normalize our results by 
((n)) as in Ref. bl) The figures below and next to the contour plot show results along the orange lines r = 3 and a = 0.6 in 
the contour plot. Solid lines are numerical results, while open circles correspond to approximation (0) taking into account the 
three pairs of singularities indicated by filled red circles in (b). (b) Singularities of the CGF in the complex plane for t — 3 
and a = 0.6. The singularities are 27r-periodic along the imaginary axis. Filled circles indicate numerical results obtained from 
Eq. (|3^), while open red and green circles show approximations ( [38[ ) and (^), respectively. Filled red circles indicate the three 
pairs of singularities entering Eq. (|l7|), as shown by open circles in (a). With time the singularities move toward the points 
indicated by filled (black) squares. 



function 



ei+a 
477(2;) 



e i+» [1 



having introduced 
r]{z) 



zl 



(34) 
(35) 



In Eq. (^) we have corrected for a missing factor of 477(2) 
in the Supporting Information in Ref. |^. This factor en- 
sures that the GF is an analytic function of z, since —77(2;) 
is the analytic continuation of rj{z) across the branch cut 
and Q{z, t) is invariant under the substitution 77 — > —77. 
Solving next for the zeros of the GF, we find 



Zk 



hi 



-,fc = l,2,.. 



(36) 



where hk is determined by the transcendental equation 

hkT = (1 + a) (/cTT — 2 arctan/ifc) , fc = 1, 2, . . . . (37) 

By graphical inspection of the transcendental equation 
we find that the solutions are real and positive for all 
k = 1,2,.... Therefore, all zeros [Eq. (p6|)] are real and 



negative, since < < 1. Remarkably, the nnaitjons 
of the zeros agree with the general statementsE3c3 by 
Abanov and Ivanov for noninteracting electrons. This 
can be understood by noting that the master equation 
[Eq. (|33|)] can also describe noninteracting spinless elec- 
trons tunneling through a resonant level for which the 
statements by Abanov and Ivanov directly apply. 

In the following, we solve the transcendental equation 
( p7| ) numerically. However, in two limiting cases, it can 
be solved analytically: For any fixed k, we can choose 
a sufficiently large r, such that hk must be small and 
arctan/ife ~ hk- We then obtain 



1 + a 

2 



k, k fixed, T large. (38) 



In the other limiting case, we fix the time t and consider a 
large k, such that hk must be large and arctan/ife ~ tt/2. 
We then find 



1 + a 



7r(fc — 1), T fixed, k large. (39) 



From the zeros of the GF we can understand the behavior 
of high-order (factorial) cumulants. Analyzing first the 
ordinary cumulants, we note that the CGF has logarith- 
mic singularities at In \zk \ + ^(21 + 1), I £ Z, correspond- 
ing to the zeros Zk of the GF. In Fig. 0b we show the 
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a) 



x10° 




b) 



5R[, 



FIG. 3. (Color online) High-order factorial cumulants in 
a function of the dimensionless time r — 2r st and the as; 
(lower plot), corresponding to the orange lines in Fig, 
to approximation (hw taking into account the three sing 
FCGF in the complex plane for r = 3 and a = 0.6. Filled 
by open circles in (a). With time the singularities move 
lie on the negative real axis, and the factorial cumulants 



0. (a) The factorial cumulant {{n''^)) 



F as 



the low- bias regime, 

mmetry a [Eq. (p2[)]. Parameters are a — 0.6 (upper plot) and r = 3 
(a). Solid lines are numerical results, while open circles correspond 
;ularities indicated by filled red circles in (b). (b) Singularities of the 
red circles indicate the three singularities entering Eq. ( |l7| ) as shown 
toward the point indicated by the filled (black) square. Singularities 
do not oscillate as functions of either time or asymmetry. 



positions of these singularities obtained from numerical 
solutions of the transcendental equation [Eq. (|37|)]. For 
comparison, we also show the limiting cases, Eqs. ( |3^ ) 
and (39), which show good agreement with the numerical 
results. In Fig. we show numerically exact results to- 
gether with the asymptotic expression [Eq. (|T7|)], taking 
into account the three pairs of complex-conjugate singu- 
larities that are closest to z = [indicated by filled red 
circles in Fig. ||(b)]. The agreement is good and the anal- 
ysis provides us with a simple interpretation of the exper- 
imental data from Ref. |[ The motion of the singularities 
in the complex plane as functions of the dimensionless 
time T and the asymmetry a cause the oscillations of the 
high-order cumulants observed in our numerical calcula- 
tions and in the experiment. 

We next turn to the corresponding factorial cumu- 
lants. Numerical results for the factorial cumulant of 
order m = 15 as a function of the dimensionless time r 
and the asymmetry a are shown in Fig. |^a and clearly 
no oscillatory behavior is observed. Again, this can be 
understood by considering the logarithmic singularities 
of the FCGF. The FCGF has logarithmic singularities at 
Zk — 1 corresponding to the (negative) zeros Zk of the GF. 
Since the singularities of the FCGF are real, the factorial 
cumulants do not oscillate as functions of either r or a 
[according to Eq. (p^]. In Fig. |b we show the positions 
of the singularities, and in Fig. 3k we show the asymp- 
totic expression [Eq. (^], taking into account the three 
singularities that are closest to z = [indicated by filled 
red circles in Fig. ||(b)]. The asymptotic expression fully 
accounts for the numerical results. 

Before closing this part, we briefly discuss the motion 
of the singularities at long times. As indicated in Figs. ||b 
and |b, the singularities of the CGF and the FCGF all 



move toward the points marked by filled black squares. 
According to Eq. (|3^), the zeros of the GF behave as 
Zk — > a^Jia?' — 1) < for long times r, since then hk — 
in Eq. (|38|). The points marked by filled black squares in 
Figs. |b and |b are thus \n[a^ / {a^ -I)] (for the CGF) and 
aV(a^- 1)-1 = l/(a2-l) (for the FCGF), respectively. 
Interestingly, the point z = a? /{a? — 1) corresponds to 
the square-root branch point of the function rj{z) in Eq. 
(^): In the long-time limit, the logarithm of the GF 
goe s as \xvQ{z,t) — j- T[rj{z) — 1]/(1 + a) according to Eq. 
(p4|), and the singularities of the CGF and the FCGF are 
then determined by the branch points of rj{z). Thus, in 
this example, all logarithmic singularities of the CGF and 
the FCGF move toward particular points in the complex 
plane, which in the long-time limit become branch point 
singularities. 



2. Finite bias 



In the general case of a finite bias, < 



< 1, 



the on-site Coulomb interaction strongly infiuences the 
charge transport statistics and we expect that the statis- 
tics will no longer be generalized binomial, with all zeros 
of the GF lying on the negative real axis. It is now a diffi- 
cult task to write down a closed-form analytic expression 
for the GF and its zeros. However, as we will see, the 
positions of the dominant zeros in the complex plane can 
be deduced from the high-order factorial cumulants. 

We first demonstrate how the zeros of the GF move 
into the complex plane due to interactions. This is il- 
lustrated in Fig. 1, where we show numerical results 
for the logarithmic singularities of the FCGF as func- 
tions of n^^\ We have solved numerically the equation 
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3 [2] 



-1 



-2 



exact result 
o asymptotics 



FIG. 4. (Color online) Motion of the dominant singularities of 
the FCGF In [Q{z + l,t)] into the complex plane. The asym- 
metry is a = 0.89 and the dimensionless time is r = 2rgi: = 4. 
Dotted arrows indicate the motion of the zeros as ng'' is in- 
creased from 0.75 to 0.85 in steps of 0.0025. Singularities start 
off on the negative real axis, but eventually they move into 
the complex plane. 



FIG. 5. (Color online) Factorial cumulants ((n™'))F as func- 
tions of the order m. Parameters are a = 0.89, r = 2Tst — 4, 
and n^*^' = 0.4, for which the dominant singularity zo — —8.2 
of the FCGF ln[Q{z-\- l,t)] is real and negative. Factorial 
cumulants have been rescaled as {{rL^})F ■ [\zo\"^/{m — 1)!]. 
Numerical results (filled squares) follow the predicted behav- 
ior for high-order factorial cumulants given by Eq. ( |l8| ) (open 
circles) . 



Q{z+l,t) = and show in Fig. ^ the singularities that are 
closest to z = 0. For low biases, n^^'' ~ 0, the dominant 
singularity is on the negative real axis, but as the bias 
increases, it is eventually reached by the second-closest 
singularity at a degeneracy point from which the pair of 
singularities splits off from the real axis and moves into 
the complex plane. At this point, the statistics is no 
longer generalized binomial due to the on-site Coulomb 
interaction. 

The change of statistics is directly reflected in the high- 
order factorial cumulants. In Fig. owe show the factorial 
cumulants ((n™)) p as functions of their order m with pa- 
rameters for which the dominant singularity is still real 
and negative, such that the high-order factorial cumu- 
lants should follow Eq. (|l|). In Fig. | we have rescaled 
the factorial cumulants as {{n"^))p-[\zo\™ / (to — 1)!], where 
Zo is the singularity closest to z = 0, and for high orders 
they clearly follow the predicted pattern with alternat- 
ing signs due to the factor (-l)(™-i) in Eq. (|lf). The 
asymptotic behavior of the rescaled high-order factorial 
cumulants (shown by open circles) agrees well with our 
numerical results (filled squares) . From the fifth factorial 
cumulant and onward, the asymptotic expression com- 
pletely accounts for our numerical results. 

In Fig. 1^ we show results for the factorial cumulants 
with parameters for which the dominant singularities 
have moved into the complex plane and the statistics 
is no longer generalized binomial. The high-order facto- 
rial cumulants are now governed by a pair of complex- 
conjugate singularities and are thus expected to follow 
Eq. (P). In Fig. | we have rescaled the factorial cu- 
mulants as {{n"^))p ■ [\zo\"^/{m — 1)!], where zg and Zq 
are the closest complex-conjugate pair of singularities. 
In this case, we see clear trigonometric oscillations as a 



function of the order m due to the factor cos (TOarg[zo]) 
with a frequency determined by arg[zo]. The asymptotic 
expression (open circles) accounts well for our numeri- 
cal results for the factorial cumulants (filled squares) al- 
ready from the third order onward. We stress that this 
oscillatory behavior of the factorial cumulants would not 
be possible in a noninteracting system for which all sin- 
gularities would lie on the negative real axis such that 
arg[zo] = TT and cos (to arg[zo]) = (—1)'™""'"'. 

For the results presented in Figs. H to |^ we found nu- 
merically the logarithmic singularities of the FCGF. In 
an actual experiment, however, the GF would typically 
not be known, and it is therefore relevant to ask if the 
position of the singularities can be deduced from the 
measured high-order factorial cumulants alone. In Ap- 
pendix ^ we describe a simple method for extracting the 
position of a dominant pair of complex-conjugate singu- 
larities from four consecutive high-order factorial cumu- 
lants. As illustrated in Fig. ^ we find excellent agree- 
ment between our numerical results and the positions 
obtained directly from the high-order cumulants (up to 
order m = 25) using this method. The figure shows re- 
sults in the full range of bias voltages from n^g^^ = to 
n|y'' ~ 1. For low bias voltages n^^'' ~ 0, the dominant 

singularity is real and negative, and as rig^'' is increased, 
two complex-conjugate singularities move into the com- 
plex plane, showing that the statistics is no longer gener- 
alized binomial. Finally, in the high-bias limit '^5'^'' = 1, 
the singularities move back onto the real axis, and the 
statistics is again generalized binomial. In the high-bias 
regime, the interaction strength U drops out of the prob- 
lem and transport takes place via the two parallel and 
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■ exact result 
o asymptotics 





15 




FIG. 6. (Color online) Factorial cumulants ((n"'))_F as func- 
tions of the order m. Parameters are a = 0.89, r = 2rst = 4, 
and n[^''' — 0.9, for which the high-order factorial cumulants 
are dominated by the pair of complex-conjugate singularities 
zq ~ —7.6 -I- 0.62 i and zq. Factorial cumulants have been 
rescaled as ((n'"))F ■ [l^o 1)!]. Numerical results (filled 
squares) follow the predicted behavior for the high-order fac- 
torial cumulants given by Eq. (|2^) (open circles). The solid 
line serves as a guide for the eye. 



FIG. 7. (Color online) Motion of the dominant singularities of 
the FCGF in the complex plane as a function of the parameter 
71^^^ The asymmetry is a = 0.89 and the dimensionless time 
is r = 2Vst = 4. Arrows indicate how the singularities move 
as n^^' is increased from to 1 in steps of 0.0025. Numerical 
solutions of the singularities are shown by filled squares, while 
open squares are the positions obtained from the high-order 
factorial cumulants (up to order m = 25) as explained in 
Appendix [BI 



uncorrelated spin channels that independently give rise 
to generalized binomial statistics. 



IV. CONCLUSIONS 

We have shown that factorial cumulants are useful for 
detecting interactions among electrons passing through 
a nanoscale device. For noninteracting electrons in a 
two-terminal conductor, the counting statistics is always 
generalizijd^ binomial, as recently found by Abanov and 
Ivanov:E£lLJ The charge transport statistics can be fac- 
torized into single-particle transfer events and the zeros 
of the GF are correspondingly real and negative. Inter- 
actions among the electrons, however, can drive the zeros 
of the GF awa^^rom the negative real axis and into the 
complex plane.E2l As we have shown, this change of the 
statistics is clearly visible in the factorial cumulants. For 
generalized binomial statistics, the factorial cumulants 
have a sign that is determined by the cumulant order 
only. In contrast, as the zeros of the GF move into the 
complex plane due to interactions, the factorial cumu- 
lants oscillate as functions of basically any parameter. 

To illustrate our findings, we have considered trans- 
port through a quantum dot weakly coupled to source 
and drain electrodes. At low bias voltages, the dot can 
only be empty or singly occupied. This corresponds to 
recent experimentsB for which we reproduce and explain 
the measured oscillations of the high-order cumulants by 
considering the zeros of the GF. In this case, the statistics 
is generalized binomial and the high-order factorial cu- 



mulants do not oscillate. As the bias voltage is increased 
an additional electron can occupy the dot. The on-site in- 
teraction now strongly affects the counting statistics, and 
the zeros of the GF move into the complex plane. As we 
have shown, this is clearly visible in the high-order facto- 
rial cumulants, which also allow us to locate the positions 
of the zeros. We expect that the motion of zeros into the 
complex plane due to interactions will be experimentally 
detectable using available measurement techniques. 

Our work leaves a number of open questions for future 
research. It would be interesting to understand physi- 
cally the exact point at which interactions cause the ze- 
ros of the GF to become complex. We have considered 
only two-terminal devices, and it might be, possible to 
generalize our findings to multilead setups. Eil 
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Appendix A: Calculation of (factorial) cumulants at 
finite times 



containing N{k + 1) elements. The equation of motion 
for \G{t))) reads 



Here we describe our method for calculating (factorial) 
cumulants at finite times without explicitly evaluating 
the (F)CGF. We begin with a master equation of the 
form given by Eq. (|28|), 



dt\g{z,t))) = M{zMz,t))), 



(Al) 



for a system with N states. Next, we substitute in 
the master equation \g{z,t))) by \g{z,t))) and M(z) by 
M(z), where \g{z,t))) = |5(e^^))) and M(z) = M(e^) 
for ordinary cumulants, and \g{z,t))) = \g{z + l,t))) and 
M(z) = M(z + 1) for factorial cumulants, respectively. 
For z = 0, the master equation then reads 



ft |g(°)(t)))=M(")|5(") (<))), 



(A2) 



having defined \g'^°'>{t))) = |ff(0,t))) and M^o) = M(0). 
Taking instead m = l,2,...,fc consecutive derivatives 
with respect to z (evaluated at z = 0), we obtain 



dtlg^'Hm 
dt\~g^'\t))) 



dtlg^'Hm 



M(i)|gW(t))) +M(0) 1.9(1) (t))), 
M(2)|g(°)(t)))+2M(i) 15^(0)) 



■3) 



(A3) 



where we have introduced the notation \g^"^\t))} = 
(9™|g(z,t)))U^o and M(") = 9;"M(z)|^^o. The (fac- 
torial) moments of order m < k are then 



(A4) 



depending on the substitutions made above. The vector 
((0| = [1, 1, . . . , 1] contains N elements equal to unity . 

We solve the system of coupled equations (A2,A3) by 
introducing the auxiliary vector 



|G(t))) = [|.g(t))),|g(i)(t))),. 



(A5) 



ft|G(0))=M|G(t))) 
where according to Eqs. (|A^,|A^) 



(A6) 



M 



M(o) J) 
M(i) M(o) J) 
M(2) 2M(i) M(o) 



(A7) 



is a matrix of di men sions N{k+1) xN{k+l). We proceed 
by solving Eq. (A6) as 



|G(t)))=eM*|G(t = 0))). 



(A8) 



Here, the initial condition as counting begins reads 



|G(i = 0))) = [|0)),0,0,...,0]^ 



(A9) 



and contains the stationary state |0)), which solves 
M(°)|0)) = 0, followed by Nk elements equal to 0. Even 
for large dimensions of M we may calculate numerically 

the matrix exponentiation eM* for a given time t and 
obtain \G{t))) via Eq. ([ASl). Having determined the (fac- 
torial) moments using Eq. (A4), the corresponding (fac- 
torial) cumulants arc obtained via the relationcll 



){F) 



m— 1 

E 

fc=i 



m — 1 
k - 1 



:("'))(F)(n'"-'= 



)iF)- 

(AlO) 



For the particular = 3 state model studied in this 
work, we could easily calculate the first m = 50 (facto- 
rial) cumulants at finite times. 



Appendix B: Determination of a pair of complex-conjugate singularities 



A pair of dominant, complex-conjugate logarithmic singularities zp and Zq can be extracted from four consecutive 
high-order factorial cumulants using methods from Refs. ^ and 52, which we repeat here for completeness. Using 
Eq. (Ell) we obtain the matrix equation 



-2 ((n- 
1 {{r: 



/Hp) 



25R(zo) 



(m 
(m 



3) 
2) 







{("" 













(Bl) 



which we solve for 3? (zq) and |zo| . Typically the accuracy of the method increases with the cumulant order m.l 
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